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Abstract 

We study an improved AdS/QCD model at finite temperature and chemical poten- 
tial. An Ansatz for the ^-function for the boundary theory allows for the derivation 
of a charged dilatonic black hole in bulk. The solution is asymptotically RN-AdS in 
the UV and AdS2 x M 3 in the IR. We discuss the thermodynamical aspects of the 
solution. The fermionic susceptibilities are shown to deviate from the free fermionic 
limits at asymptotic temperatures despite the asymptotically free nature of the gauge 
coupling at the boundary. The Polyakov line, the temporal and spatial string tensions 
dependence on both temperature and chemical potential are also discussed. 



1 Introduction 



Non-perturbative QCD at finite chemical potential is still elusive. First principle formula- 
tions such as the lattice suffer from the sign problem [1] . Most non-perturbative formulations 
based on semiclassics such as the instanton or dyon formulations require further insights on 
the role of the fermionic zero modes at finite chemical potential. Some insights on the role 
of the chemical potential at finite temperature can be gleaned from strong coupling lattice 
QCD [2] or models [3]. Most of these models for light 2-flavor QCD suggest a second order 
transition at small chemical potential and finite temperature, and a first order transition at 
higher chemical potential. Noteworthy is the occurence of a tricritical point which appears 
to be sensitive to the nature of the confining forces. 

Non-perturbative QCD with a large number of colors at finite chemical potential is likely a 
crystal of confined baryons. The crystal binding energy can be parametrically small, causing 
it to melt under quantum fluctuations and/or temperature. The result is a strongly coupled 
baryonic liquid. The original Skyrme model supports this descriptive although it suffers 
from the inherent shortcomings of the higher order chiral terms at high density [4]. The 
chiral holographic approaches to QCD support the crystal structure in the large number of 
colors limit without the shortcomings of the Skyrme model [5]. The crystal is found to melt 
at relatively small temperatures by the Lindemann criterion, resulting into a holographic 
liquid of instantons at low density and dyons at higher densities. Therefore, it is reasonable 
to think that the cold and dense inhomogeneous phase gives way to a homogeneous phase 
under the effects of temperature. 

The purpose of this work is to address issues in relation to the dense and hot homogeneous 
baryonic phase using the (bottom-up) holographic approach. Specifically, we will use a 
variant of the improved (bottom-up) holographic approach to finite temperature Yang-Mills 
in the large number of colors limit put forward by Gursoy et al. [6, 7, 8] with the addition 
of a fermion chemical potential to account for baryon number. In [6, 7] a potential for 
the dilaton field is constructed to reproduce some key features of Yang-Mills theory in 4 
dimensions, namely heavy quark confinement in the infrared and asymptotic freedom in the 
ultraviolet. Recently, a variant of this construction was suggested in [8] whereby the dilaton 
field is directly tied to the running of the gauge coupling in the Yang-Mills theory. This 
construction is more transparent physically as it directly ties the holographic direction to 
the running of the gauge coupling constant. It is also less numerically intensive. We follow [8] 
and add a U(l) charge field in bulk to account for the effects of a finite chemical potential 
at the boundary. 
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In section 2 we discuss the bottom-up improved holographic model in 5 dimensions with 
bulk gravity coupled to a dilaton and the U(l) charged field. The dilaton dynamics follows 
from the running of the gauge coupling. The dilaton potential is fixed by the equations of 
motion. Explicit solutions are constructed in section 3. In section 4 we discuss the small 
and big charged black hole solutions, and their ensuing bulk thermodynamics. In section 5 
we analyze the fermionic susceptibilities and compare them to the ideal gas limit. In section 
6, the Polyakov line is evaluated at finite temperature and chemical potential. In section 7 
we discuss the role of the chemical potential on the spatial and temporal string tension. Our 
conclusions are in section 8. The sensitivity of the model to the dilaton-gauge-field coupling 
is discussed in Appendix A. In Appendix B we derive (quark) susceptibilities for a warped 
holographic model without a dilaton. 



2 The Model 

In the Einstein frame the 5-dimensional Einstein-Maxwell action in the background of a 
scalar or dilaton field 0, is given by [6, 7] 

S = d 5 xV=g(R - \d A <t>d A <t> + V{4>) - ^F AB F AB ) - 2 J d'xVhK^j (1) 

with the U(l) gauge field tensor Fab = <9[a^4_b] and L the radius of the AdS-space in 
the conformal limit. The Gibbons-Hawking boundary term gives no contribution to the 
equations of motion but is crucial for evaluating the on-shell action [9]. The coupling to the 
£7(1) charge is a generalization of [10], where we consider the parameter c in the exponential 
gauge coupling as a free parameter. A recent analysis in this direction is found in [11]. For a 
non-constant scalar field <fi = <f>(z) the Bianchi identity V a Gab = (with Gab the Einstein 
tensor) ensures that the Einstein equations imply the equation of motion for the scalar field. 
With the following Ansatz for the metric 

ds 2 = b(z) 2 (-f(z)dt 2 + dx 2 + j^ dz ^ (2) 

we follow [6, 7, 8] and assume that the /3- function for the gauge field theory on the boundary 
at z = is given by 

m = b^ = -a>a« (3) 
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with the running t'Hooft coupling X(z) = ~ g\ M N c . Note that (3) does not follow from 
varying the action (1). q > 1 ensures confinement in the IR and the values q = 10/3 and 
Po = 488.8 have been shown in [12] to reproduce the dilaton potential in [13] to lowest order 
in A. We will use these values in our numerical analysis. For a static, charged black hole 
solution the equation of motion for the scalar potential A (z) is given by 



b( Z y 5 d z (b(z)e c ^A' (z)) oc S( 



Z ^charge ) 
p p -c(p(z) 

with A' (z) = -j^A (z) and the electric charge e located at some position z c harge behind the 
horizon of the black hole. We can now solve the Einstein equations 

b' 2 b" 4 

4 - 4 - 1^ - < 5 > 

§x + ^-^ = ° (6) 
9 6'/' jf , / &v n m 

together with (3) and (4). Note that the factor f(z) drops out from the spatial components 
of the energy- momentum tensor for the gauge field. 



3 Solutions 

The solution for the electrostatic potential is given by 

e f z e~ c ^ x ^ 

Mx) = --^ dx — r+ll . (8) 

The integration constant \x is interpreted as the chemical potential since A (z — > 0) = ji [14]. 
The requirement for the electrostatic potential to vanish at the horizon of the black hole, 
see (23), gives 

Mz H ) = 0*/. = ^ d X —. (9) 

The warping factor b(z) is not influenced by the presence of the charge and the solution to 
(3) is given by 



0=h iT<FT5AAF*- (10) 
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Since 

W(X) = ^ = W(0)e% (11) 



z = ^ dx e~ x -^ , (12) 



and a = {w=qj) ' ( 5 ) y ields 

b JQ 

where the constant l/W(ti) = L is fixed by the boundary value of the potential V(z) ) see 
(22) and [8]. (10) and (12) imply that neither the warping b(z), nor the coupling A depend 
on the temperature or chemical potential. (12) introduces a scale A = ^ for our model and 
implies an upper bound on the radial coordinate 



£ 



oc 



z < — dxe x 4^ — Zw . (13) 
bo Jo 

A static wall in the IR is generated which leads to an area law for the Polyakov loop as we 
will detail below. 

Expanding the integral in (12) for small z (large Q) 

/*00 < ^~ > 1 

J dxe-°-* = £_(-£) B r(l-n,Q) (14) 



o 

exp 



-0- — 
W 4Q 



+ (15) 



with the incomplete Gamma-function T(l — n, Q) we obtain b(z — > 0) = - and (f>(z — > 
0) oc — log (— log Az), which reproduces a logarithmic running coupling for the model on the 
boundary. 

Changing variables from z to Q in (6) reads 

%f(Q) + (~ + 4) *,/(«> = (fa-W^ e3 - e0QA (16) 

The two asymptotic conditions on f(z), f(z = 0) = f(Q = 00) = 1 and }'{Qh) = at the 
horizon Qh = Q(zh), fix the two integration constants and we obtain 

/(g) = 1 + Cl i(Q) + 5^% [ Q dx e-^" 4 * f dy e'^y^ (17) 

JlV-ffJ Joe Joe 
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with 5 = — ^- and 

_1 _ 5 ^fX f QH dx e-^- 4 * P dy e-^- 2y y^ 
j(Q) = / e-£ _2x a;^5 (19) 

J oo 

i(Q) = / dxe~^- 4x . (20) 

We see that our solutions for f(z) reduces to the one obtained in [8] for /j, — 0. With a 
vanishing coupling of the dilaton to the gauge field, i.e. c = 0, the metric in the Einstein 
frame (2) approaches the RN-AdS-metric for z — > with 

/(^O^l-^Atf + Jg^)'. (21) 

We can now solve (7) for the scalar potential 

V(Q) = W(0fe^ {l# Q + (5 - ^)d Q + 12 - ^)f(Q) (22) 

As is temperature independent this implies V(Q) = V(<f>,T, [/,). The dependence on tem- 
perature and chemical potential follows from f(Q). The normalization V(z — > 0) = V(Q — > 
oo) = § gives ^ = L. 

4 Thermodynamics 

The Hawking temperature of the black hole, 

T = ±f(z H ) , (23) 



is given by 



" 4vr i(Q H ) V ' M/ ' A) j(^) 2 1 



with z H , Q H = Q(zh) the position of the horizon and 

"0 



k(Q) = / dxe~^- 4x / dye-^~ 2y y^. (25) 

7 oo •/ a; 
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Figure 1: Left: The Hawking temperature at chemical potentials /i/A = (solid), 1 (dashed), 
3.5 (dotted) as a function of the scaled horizon z#A. Right: The chemical potential at T = 0. 



At vanishing chemical potential the black hole solution has a minimal non-zero tempera- 
ture. The two regions T\z H ) > (T'(Q H ) < 0) and T'(z H ) < (T'(Q H ) > 0) correspond to 
a small and a big black hole branch. Since the solution for \i = shows a minimum temper- 
ature, a certain density (/i = 2.8A) is needed to obtain a solution with zero temperature. 
Here, we will focus on the case of vanishing dilaton-gauge-field coupling, i.e. c = 0, and 
consider the case c = 4 separately in Appendix A. For T — > and /i > /xq, the small black 
hole vanishes leaving only the big black hole solution, see Fig. 1, and the metric in the IR 
becomes AdS 2 x IR 3 , [15, 16]. With (z - z H ) = t/C, * = 2/{lf"{z H ))T and 7 ->■ 0, C, r finite, 
(2) reduces to 



One crucial step in the above analysis is that the temperature as a function of the horizon 
does not have a minimum as T — > 0, i.e. /" {zh{T = 0,/i)) 7^ 0. In our case this translates 
to the disappearance of the small black hole as T — > 0. 

Note that regularizing the action in (1) by subtracting a 'vacuum' (thermal gas) solution 
with functions <po,bo, fo = 1 does not necessarily yield the grand potential since the scalar 
potential (22) has additional temperature and chemical potential dependence thereby up- 
setting the Gibbs relations. This notwithstanding, we can still define the entropy density 
s(T,n) carried by the black hole through its area as 



(c = 0) 



ds 




(26) 



s(T, fj,) 



- "0 r 3Q H 



(27) 
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The small (big) black hole branches have negative (positive) specific heat, c v oc Tds/dT, and 
are unstable (stable). 



r 3 




1.5 2.0 2.5 To 



Figure 2: The scaled entropy density s/T 3 at chemical potentials fi/A = (solid), 1 (dashed), 
1.5 (dotted) with To = 170MeV and A = 321MeV. Both the contributions from the small 
and big black hole are shown. 



The scaled entropy density s/T 3 reaches its asymptotic value around T ~ 2T and de- 
velops a peak with increasing \i. A low temperature and high density expansion yields 

bf.fl u 3 7tu 2 T n 2 3VS uT 2 anVS l±T 2 n/rT1 . A , . 

4G 5 A 3 4 A 3 Sv^ A 3 16^ A 3 log ^ '^y 

(28) reduces to the RN-AdS entropy [17] up to logarithmic corrections. Much like the RN- 
AdS black hole, our black hole solution carries a finite entropy at zero temperature and a 
linear specific heat, see Fig. 4. The effects of the running coupling through the scalar field 
causes only logarithmic corrections much like the high temperature case. 

The pressure at vanishing chemical potential is obtained by integrating the entropy den- 
sity over the small and big black hole branch 

p(T,fi = 0) = dQ - s(Q(T,fi = 0)). (29) 

The scale A = ^ is fixed by requiring that the pressure at /i — vanishes at a critical 
temperature T = 170MeV or A = 321MeV. Comparison with an ideal gluon gas gives 

_ 16iV e 2 A 3 r R i 
G 5 45vr L°J • 

The pressure of the small black hole is negative, indicating its instability. Fig. 3 shows the 
big black hole contribution to the pressure. For T/T < 1 the 'vacuum' solution dominates. 
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The integration constant in (29) is fixed by analyzing the high temperature thermodynamics 
of the small black hole, [9]. Since the small black hole vanishes for T = 0, we fix the 
integration constant by setting the zero temperature pressure to zero at fi = fi min 
n(T = 0,Q H = 0) =2.8A, 



T=0 



p(T = 0,fi) 



rCr=M 



dQ[^)n(Q{T = 0,ti)), 



(30) 



with the charge density n given by (31). Fig. 3 shows that the black hole solution is 
unstable against the vacuum solution below /i . In contrast, the RN-AdS black hole at zero 
temperature [18], dominates the action for all chemical potentials down to \x = 0. 




P(0. H) 



T 
To 




Figure 3: Left: Big black hole pressure at pi = (N c = 3) and pressure of an ideal gluon gas 
Pideai gas /T A = 2 {N 2 C - 1) g (dashed). Right: Pressure at T = (N c = N f 
ideal relativistic Fermi gas of Nf massless quark flavors Pideai gas/ ^ 



N f 



For fixed chemical potential the charge density n(T, fi) is proportional to the electric 
charge e and using (9) we write 

a aL 3 fx — a5 _1 / 3 /i 

with the proportionality constant a and 

m (Q) — J dxExp ^— 2xj x~ . (32) 



8 



For /i>Twe obtain 



2a (la 3 , 
n(/i»T) = - -£r(l + 



+ 



/1 2 T 7T 



aiv(3 + \/6) 



A \6A3^ 4log-fe' A3 ^ 241og^ 



) + 0(T//i) . (33) 



As in the RN-AdS, the leading order low temperature correction is linear in T. The deviation 
from the pure RN-AdS shows up as logarithmic corrections. 
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Figure 4: Left: Scaled charge density n/T 3 for the big black hole at densities fi/A = 
0.5 (solid), 1 (dashed), 1.5 (dotted) with a = ^2 A 4 . Right: s/n 3 (solid), nf/j, 3 (dashed), 
e//x 4 (dotted) at T = 0. 



5 Susceptibilities 

In QCD at low density and high temperature the pressure can be expanded as 

P(T, /i) 



rp4 



5>( T )(f) n (34) 



n=0 



with the (flavor symmetric) quark susceptibilities 

c (T) - LJ^-VSElA (35) 

Various hadronic susceptibilities in the transition region around To = 170MeV at vanishing 
quark chemical potential show distinct characteristics. In particular, lattice studies [19] 
and expectations from PNJL models [20] both confirm that c 4 , c 6 show distinct peaks at a 
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critical temperature T c . Asymptotics of the susceptibilities come close to the ideal gas values 
at temperatures T ~ 2T C . These charge density fluctuations are obtained as derivatives with 
respect to the density on the grand canonical partition function. 

The quark number susceptibility in hard and soft wall AdS/QCD models were studied in 
[21, 22]. In these models the first non- vanishing coefficient in the expansion (34) show a jump 
at a critical temperature due to a Hawking-Page transition. In the improved model under 
consideration with running gauge coupling, we can explicitly assess the first few moments in 
(34). We first note that the odd moments vanish since m(Q H ) in (31) receives corrections 
at high T (high Q H ) of the form /j, 2 /T 2 . From (31) we obtain 



Note that eqn. (23) implies 



c 2 (T) 



1 d 

-ad' 1 ' 3 



n{T,n) 
1 



/i=0 



2T 2 A 2 m(Q H ) 



At=0 



OQ 



H 



dfx 



= , 



(36) 
(37) 

(38) 



d 2 Q H 
dfj, 2 

d 4 Q H 
d/i 4 



n=o 



We obtain for c 4 (T), c 6 (T) 



Exp (~3Q H ) 



k{Q 



H 



2nAT 



3i(Q H ) + Exp{-AQ H - ^-) 



-3 



d 2 Q H 



d 2 T 



9Qh 

or 



c 4 (T) = 



1 d 3 

-a5- l ' ? > 

24A 2 
T 2 d 5 



n(T, fj) 



d 2 Q H \ d 



6! (9/i 5 
6!A 2 



dfj, 2 
n(T, fj) 



dQ H m(Q H ) 



li=0 



A 



4 Qh\ d 



JdQ H 



U(Qh)) 2 ,=o ' 



fi=0 



d 2 Q H \ 2 d 



(39) 
(40) 



dfi 2 ) dQ\) m(Q H ) 



fi=0 



(41) 
(42) 
(43) 
(44) 



As can be seen from a high temperature expansion of the integrals in (37), (44) (Q H — > oo) 
C4 and cq do not converge to a finite asymptotic high temperature value unless c = 0. Thus, 
a meaningful comparison to QCD demands a vanishing coupling of the dilaton to the U(l) 
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charge. With c = the high temperature asymptotics are given by 





As a result of the non-trivial warping b(z), the susceptibilities are non-constant and a 
comparison with recent lattice data [19] shows that the susceptibilities C2, C4, Cq obtained in 
this model have the correct shape around the critical temperature To. C2 and c§ approach 
its asymptotic value from below, while C4 shows a distinct peak at the critical temperature. 
All high temperature asymptotics of the susceptibilities are strictly positive. For the model 
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proposed in [23], the susceptibilities vanish as explained in Appendix B. 



In comparison, for an ideal gas of massless quarks and anti-quarks the susceptibilities read 

jdeal gas = AUV^ jdeal gas = Ag, ? jdeal gas = g ^ thg jdeal gas/geal gas = l/^7T 2 ), 

Cq/c 4 = 0. In our model: c 4 /c 2 = l/(187r 2 ), c 6 /c 4 = l/(307r 2 ). While the scale A was fixed 
by the gluonic part of the pressure, we fix the constant a by comparing to the fermionic part 
of the pressure of an ideal gas. For N c = 3, Nf = 3 we fit the high temperature asymptotic 
of c 2 to the ideal gas value and obtain a = ^A 4 . The susceptibilities in the holographic 
model considered in [17] with a non-confining RN-AdS black hole are constant and given 
by c* N ~ AdS = iV c 2 7 2 /8, cf N - AdS = gf, cf N ~ AdS = -gj, where the flavor dependence is 
embedded in the parameter 7 2 oc Nf/N c . 




It is now straightforward to evaluate the density contribution to the pressure, Ap(T, \i ^ 
0)/T 4 = Er= 2 Cn(T) (f ) n ~ c 2 (T) (§) 2 + c 4 (T) + c 6 (T) (^) 6 , and the energy density 
e(T,fj) — Ts — p + /in at high temperatures and small charge density. Fig. 7 shows our 
results for (e — 3p)/T 4 at different densities. 

6 Polyakov line 

The expectation value of the Polyakov line can be schematically written as [24, 25] 

(L(T,»)) = Y,*>nExp(-S n ) , (48) 

n 
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with weights w n associated with a renormalized area S n . We will approximate the area using 
the Nambu-Goto action describing a fundamental string stretched between the horizon and 
the boundary at z — 0. While on-shell quantities such as the susceptibilities do not depend 
on the frame, quantities such as the Polyakov line may depend on it. The appropriate 
background for the string is given by the metric in the string frame 

9ab(*) = e^ z) g AB {z) , (49) 

where qab is given in eqn. (2). For a static configuration with the parametrization £i = t, 
£ 2 = z the action is given by 

Sng = ^7 J d 2 ^g\ B d M X^d N X^ (50) 

= J dzb(z) 2 \(z)sVTTJ(zWW- (51) 

The first integral gives the equation of motion for x(z) 

d ( b(z) 2 X(z)lf(z)x'(z 



dz \ y/TTJWW 



0. (52) 



It is easy to check that the dominant contribution to the integral in (51) comes from the 
solution x = const., and we will neglect other solutions. Since our metric is asymptotically 
RN-AdS, implying that b(z — > 0) = ^, the action defined in (51) is divergent. We renormalize 
the action by subtracting the vacuum contribution stretching from the boundary (z = 0) 
upto the wall at z w . The Polyakov line in the saddle-point approximation (48) reads 

(L(7»> ~ w Exp(-S r N e Z) (53) 

(1 f ZH 1 f Zw 4\ 

dzb ^ 2X ^~ 3 -^fJ o dzb(z) 2 \(z)3) (54) 

= w » Ex v{^f dzb(z) 2 \(z)iy (55) 
At vanishing chemical potential, the 'low' temperature limit of the Polyakov line is given by 



4 



(HT -> r m „,o)> « e xp EXP ( -x^-j) (x^-y) *'~" } (56) 
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and, thus, shows little deviation from its high temperature asymptotic value. Fig. 8 shows 
the behavior of the approximate order parameter of the phase transition. The constant 
b\ja' is fixed by comparing the spatial string tension with lattice data, see next section. At 
fixed \i the expectation value of the Polyakov line jumps rapidly from a minimum value at 
a temperature T = T min and approaches a constant value in the high temperature region. 
A direct comparison with lattice data is unfortunately not possible given the subtraction 
dependence inherent in the definition of the Polyakov line both on the lattice and in our 
case. This point is generally overlooked in most analyses. 



T 




Figure 8: Left: Polyakov line with /i/A = (solid), 1 (dashed), 1.5 (dotted); (w = 1). 
Right: Big black hole contribution to the spatial string tension at densities /i/A = 
(solid), 1.5 (dashed), 2.5 (dotted). 



7 String tensions 

The analyses of the spatial and temporal (effective) string tension, a s and a at finite tem- 
perature and density, are identical to those carried in [12] at finite temperature. Indeed, a 
rerun of their analysis shows that for the spatial string tension a s , the density dependence 
is entirely encoded in the position of the horizon zh(T, /i) with 



2na 



Z = Z H 



2na 



-b(z) 2 X(z)s 



Z = Z H 



(57) 



where b s (z) is the warping in the string frame. A comparison with lattice data [26] at T = To, 
/j, — fixes 4% = 3 * 10" 4 MeV _1 . For chemical potentials up to ji — 2.5A, the spatial string 
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tension shows a weak dependence on the chemical potential in this improved holographic 
model. 

In the temporal Wilson loop, the string tying a heavy quark and antiquark in the bi- 
fundamental representations extends in the holographic direction in bulk. The further the 
spatial separation L, the deeper the string extends in bulk z. Specifically [12] 



r i 

L(z*) = 2 dz (58) 

with z* the maximum holographic depth for the pending string. In the black hole background 
with / 7^ 1, L(z*) is finite for any temperature and chemical potential as shown in Fig. 9. 

LA 

, LA 




Figure 9: Length separating two heavy quarks. Left: /i = with T = T (solid) and 
T = 1.5T (dashed). Right: T = T , n = /i = 2.8A (solid) and T = 0,/i = 4A (dashed). 



In the black hole background two heavy-quarks still tie up to distances of order L max 
which is the location of the maximum in Fig. 9. In other words, for < L < L max we expect 
a linear-like free energy between two heavy quarks, while for L > L max the heavy quarks 
"screen" through the large entropy of the string (at zero density) or through fundamental 
matter (at finite density) causing the free energy to plateau at twice the screening masses. 
Around the phase transition points, T = T and fi = 0, L max = 1.24A -1 ps 0.76 fm, while 
for T = and /i — /iq — 2.8A, L max = 0.10A -1 ps 0.06 fm. The larger the density, the larger 
the screening. 
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8 Conclusions 



Finite density QCD both at zero and finite temperature is still a challenging problem from 
first principles. In the double limit of a large number of colors and large t'Hooft coupling 
the holographic approach offers a non-perturbative tool for investigating QCD-like gauge 
theories at finite temperature and density. The current analysis provides a step in that 
direction whereby finite density effects are incorporated in an improved AdS/QCD model 
with a bulk U(l) charge. Without a dilaton-gauge coupling, our improved AdS/QCD metric 
at finite (charge) density is asymptotically RN-AdS on the boundary with a non-trivial 
dilaton profile interpreted as a running coupling. The effects of the density do not alter the 
warping factor of the underlying gravitational metric. 

At zero density but for temperatures larger than a minimum temperature, the gravita- 
tional equations yield a pair of black holes of different holographic sizes. A small black hole 
that is unstable thermodynamically, and a large black hole that is stable. The occurence of 
a minimum temperature in the improved model reflects on the (first order) transition from a 
'vacuum' (thermal gas) to a black hole solution. The unstable solution is found to disappear 
at large densities. For zero temperature, the small black hole solution disappears while the 
large black hole solution requires a minimum chemical potential /iq = 2.8A. As a result, the 
metric becomes AdS 2 x M 3 in the infrared. Evaluating the pressure at finite /x and T = 
shows that the 'vacuum' becomes unstable against black hole solution at a critical \i = /i 
indicating a first order Hawking-Page transition. 

The (quark) susceptibilities C2, C4, cq show rapid variations across the critical temperature 
T in the improved holographic model, in contrast to the RN-AdS model where they are 
found to be constant. The variations are overall consistent with the ones reported by lattice 
simulations. However, their asymptotic ratios are off compared to the free fermion (quark) 
limits, despite the fact that in the improved holographic model the gauge coupling on the 
boundary runs weak at high temperature. 

To further clarify the nature of the charged black hole solutions, we have analyzed the 
subtracted Polyakov line and found that it jumps at both finite temperature and/or finite 
density. The jump is suggestive of screening which is expected to be larger at finite charge 
density. Indeed, our analysis of the temporal Wilson loop confims that two heavy quarks 
detach more easily at higher density. Spatial Wilson loops are midly affected by the charge 
screening in the black-hole background. 

The bottom-up model with running coupling constant appears to capture some essentials 
of a first order transition from a vacuum with matter characterized by / = 1 and a charged 
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black-hole characterized by / ^ 1. This first order geometrical transition captured by the 
holographic gravity equations appears to encode some of the features expected from a first 
order transition in QCD with a large number of colors in the homogeneous regime. The 
latter is expected to be dominant for a broad range of temperatures and densities. However, 
it has two major shortcomings: 1/ The (quark) susceptibilities asymptote the wrong values 
despite the fact that the bulk thermodynamics (pressure, entropy and energy densities) 
can be adjusted to asymptote the correct values. 2/ The very low temperature phase is 
likely inhomogeneous. While the latter issue is readily overcome [5], the former issue is more 
problematic at higher temperature. It goes to the heart of high fermionic charge fluctuations 
in dense matter and therfore the very interpretation of our charged black hole solutions. 
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A Non-trivial dilaton-gauge coupling: c = 4 



In this Appendix, we show that a non-trivial dilaton-gauge coupling affects considerably the 
formation of the large black hole solution in the model we discussed above. 
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Figure 10: Left: The Hawking temperature at chemical potentials /i/A = 
(solid), 1 (dashed), 3.5 (dotted) as a function of the scaled horizon z H A. Right: The 
density at T = 0. 



Fig. 10 shows, that the big black hole branch vanishes completely for densities larger 
than fi leaving only the unstable solution with an absolute maximum temperature in the 
IR region (zh —> z w ). In this region, the size of the black hole at a given temperature is 
insensitive to changes in the charge density. Discarding the unstable small black hole, the 
range of validity for this setup is dictated by T > T min , ji < /i Q . 

s 



p(0, fi) 




20 40 60 80 To [ 

Figure 11: The scaled entropy density s/T 3 at chemical potentials fi/A = 
(solid), 1 (dashed), 1.5 (dotted) and the pressure at T = 0. 



The scaled entropy density, s/T 3 , peaks at finite /i around T = 5T and the high temper- 
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ature asymptotic behavior is much slower as in the case of vanishing dilaton-gauge coupling, 
Fig. 11. The susceptibilities C4, c$ do not asymptote for large temperatures. The pressure at 
T = is negative for small /i/A and asymptotes to zero for /i/A > 10. Thus, the 'vacuum' 
solution is stable against the black hole solution at T = for all ranges of the chemical 
potential and no phase transition expected. 



T 
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Figure 12: Left: Polyakov line with /i/A = (solid), 1 (dashed), 1.5 (dotted). 
Right: Big black hole contribution to the spatial string tension at densities /i/A = 
(solid), 1.5 (dashed), 2.5 (dotted). (w = 1) 



The Polyakov line shows a peak around T = T that is enhanced with increasing charge 
density. Fig. 12 shows the results for the Polyakov line and the spatial string tension with 
non-trivial dilaton-gauge coupling. 

B Warped Model 

We discuss the asymptotics of the susceptibilites as obtained for the warped metric proposed 
in [23] with a charged black hole. At finite temperature and density the metric reads 

rfs 2 = ^ e o.45GoW (-f( z ) d1 ? + dx 2 + ^--dz 2 \ (59) 

z 2 V /(*) / 

and we assume that the relation between the horizon zh and the temperature, density is 
given by the RN-AdS result [18] 



:„ =2 I kT+ \ TT^-lyV ) ■ (CO) 
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The solution (9) is generic and we obtain 

„ — _£_/ ^ r ^ (3 -0.225GeV 2 x 2 e 0.225GeV 2 2 ? J / fil \ 

'""Wo to £ e " £4 . 45GeV 2 e • ^) 

The charge density is proportional to the charge and using (61) yields 

n(T,fi)ocfie°- 225GcY2z «. (62) 

Unlike the improved holographic model, the corresponding susceptibilities vanish as power 
laws at high temperature, 

1 0.225GcV 2 

c 2 oc — e^^- -> (63) 
1 z 

1 0.225GeV 2 

c 4 ex — e^^- -> (64) 

(1 1 \ 0.225GcV 2 

#— + #— Je^s-^0 (65) 

This warped and charged black hole cannot be used to model up the fermionic fluctuations 
in screened QCD. 
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